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Abstract. We consider the Casimir force acting on a d-dimensional rectan- 
gular piston due to massless scalar field with periodic, Dirichlet and Neumann 
boundary conditions and electromagnetic field with perfect electric conductor 
and perfect magnetic conductor boundary conditions. The Casimir energy in a 
rectangular cavity is derived using cut-off method. It is shown that the diver- 
gent part of the Casimir energy does not contribute to the Casimir force acting 
on the piston, thus render an unambiguously defined Casimir force acting on 
the piston. At any temperature, it is found that the Casimir force acting on 
the piston increases from — oo to when the separation a between the piston 
and the opposite wall increases from to oo. This implies that the Casimir 
force is always an attractive force pulling the piston towards the closer wall, 
and the magnitude of the force gets larger as the separation a gets smaller. 
Explicit exact expressions for the Casimir force for small and large plate sep- 
arations and for low and high temperatures are computed. The limits of the 
Casimir force acting on the piston when some pairs of transversal plates are 
large are also derived. An interesting result regarding the influence of temper- 
ature is that in contrast to the conventional result that the leading term of the 
Casimir force acting on a wall of a rectangular cavity at high temperature is 
the Stefan-Boltzmann (or black body radiation) term which is of order T d+1 , 
it is found that the contributions of this term from the two regions separating 
the piston cancel with each other in the case of piston. The high tempera- 
ture leading order term of the Casimir force acting on the piston is of order 
T, which shows that the Casimir force has a nontrivial classical h — > limit. 
Explicit formulas for the classical limit are computed. 



1. Introduction 

In 1948, Casimir predicted the existence of an attractive force between two per- 
fectly conducting parallel plates, which is due to the zero-point fluctuation of elec- 
tromagnetic field between the plates pQ. Nowadays, Casimir effect is generally 
referred to similar effect due to any quantum fields. Since 1948, thousands of pa- 
pers due to Casimir effect have appeared in the literature, and it has gained more 
and more attraction from physicists and even engineers. Many experiments have 
been designed to verify the existence of the Casimir force. For example, Mohideen 
et al [2] used atomic force microscope to confirm (within a few percent error) the 
existence of Casimir force for plate-sphere separation between lOOnm to 900 nm. 
At these length scales, Casimir force becomes non-negligible and therefore people 
working in nanoscience and nanotechnology starts to show interest in this effect. 



Key words and phrases. Field theories in higher dimensions, finite temperature Casimir effect, 
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Research has been done on how to use Casimir force to drive a nanodevice, as well 
as to eliminate unpleasant effect due to Casimir force, such as adhesion or stiction. 

In the conventional calculation of Casimir energy, one sums over different modes 
of zero-point energies in the presence of boundaries: 



However, this sum is divergent and regularization methods are used to obtain a 
finite result. One way of regularization is by subtracting the vacuum energy due 
to the background infinite space. However, this procedure does not always lead 
to a finite quantity, due to the possible hypersurface divergence. As has been 
pointed out by Deutsch and Candelas [5] and Baacke and Krusemann 4 , it is 
impossible to compute the Casimir energy with boundary conditions to approximate 
the interaction of the vacuum fluctuations with the boundary material. This has 
been further discussed by Graham, Jaffe et al in a series of papers El El [5J . They 
argued that the surface divergences cannot be removed by any renormalization 
of the physical parameters of the theory. In 2004, Cavalcanti [9] showed that 
despite the presence of the divergence in the Casimir energy, it is possible to obtain 
unambiguous finite Casimir force for a special geometric setup known as piston. 
He showed that for a two dimensional massless scalar field theory in a rectangular 
piston, the surface divergent terms of the Casimir force on the piston due to the 
two regions divided by the piston cancel each other and the resulting Casimir force 
acting on the piston is always attractive. Since the original work of Cavalcanti, 
piston geometry and its variants have attracted considerable interest. In [TUl [Tl] , it 
was shown that for electromagnetic field with perfect electric conductor conditions 
inside three dimensional rectangular piston, the surface divergent terms are also 
canceled and the resulting Casimir force acting on the piston is attractive. The 
Casimir piston for three dimensional electromagnetic field with prefect conductor 
conditions are studied further in [T21 [T31 [H] , where pistons with arbitrary cross 
sections are considered. In p~5], the rectangular piston for massless scalar field 
with Dirichlct boundary conditions in three dimensions was discussed in detail. 
In [16], it was proved that for massless scalar field with Dirichlet and Neumann 
boundary conditions, the hypersurface divergent terms of the Casimir force due to 
the two regions separating a <i-dimensional rectangular piston always cancel each 
other, and the Casimir force is always attractive. The Casimir force on rectangular 
piston due to electromagnetic field with perfect magnetic conductor conditions are 
computed and discussed by Edery and Marachevsky [17] . In another recent work 
[15] . the Casimir force for massless scalar field with Dirichlet boundary conditions 
on a rectangular piston in the spacetime with extra compactified dimensions was 
discussed. In [19] , it was shown that the Casimir force between two (nonmagnetic) 
dielectric bodies which are related by reflection is always attractive. This result was 
generalized by Bachas [20] who showed that reflection positivity implies that the 
force between any mirror pair of charge-conjugate probes of the quantum vacuum 
is attractive. The attractive nature of the Casimir force will create undesirable 
effects such as the collapse of a nano device - an effect known as stiction [55] . 
Therefore, it becomes desirable to search for circumstances where the Casimir force 
can be made less attractive, or even repulsive. In [23], Barton showed that for a 
thin piston with weakly reflecting dielectrics, the Casimir force at small separations 
is attractive, but turn to repulsive as the separation increases. In [24], it was shown 
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that in one, two or three dimensions, if one surface assumes Dirichlct boundary 
condition and the other assumes Neumann boundary condition, then the Casimir 
force on a rectangular piston is repulsive. Another scenario that leads to repulsive 
Casimir force was discussed in |25j . We would like to remark that a scenario similar 
to piston geometry has been considered by Reuter and Dittrich in 1985 [26] in 
the context as a mechanism to regularize the Casimir energy. Recently, a similar 
mechanism is used to compute the Dirichlet Casimir effect for </> 4 theory in (3+1) 
dimensions [2"T] . 

In this paper, we study the Casimir effect of massless scalar field and electro- 
magnetic field in a d-dimensional rectangular piston at finite temperature. Differ- 
ent boundary conditions are considered including periodic, Dirichlct and Neumann 
boundary conditions for massless scalar field and perfect electric conductor (PEC) 
and perfect magnetic conductor (PMC) boundary conditions for electromagnetic 
field. In this paper, we propose an alternative method to calculate the cut-off de- 
pendent Casimir energy. We notice that the divergent part of the Casimir energy 
is linear in each of the variables Li, . . . , Ld — the lengths of the rectangular cavity. 
This immediately implies that the divergent terms do not contribute to the Casimir 
force acting on the piston. Therefore the Casimir force acting on the piston is finite 
and unambiguous without any renormalization. The thermal effect of the Casimir 
force has so far only been considered for electromagnetic field with PEC boundary 
conditions in three dimensions [TTI [TBI IT4] , However, [TTJ does not give a correct 
high temperature behavior of the Casimir force. In this paper, we derive for d- 
dimensional piston geometry explicit and exact formulas of the Casimir force which 
can be used to study the high temperature and low temperature behaviors of the 
Casimir force. We show that in the high temperature regime, the black body radi- 
ation contributions to the Casimir force from the two regions separating the piston 
cancel with each other. The leading order term of the Casimir force is of order T, 
which shows that the Casimir force has a classical limit. 

In this paper, we choose the units where h = c = k = 1. 



2. Casimir energy by exponential cut-off method 

Consider a <i-dimensional rectangular piston, which is a rectangular cavity [0, L\] x 
[0, L2] x ... [0, Ld] separated by a hyperplane x\ = a (the piston) to two regions: 
the region I of dimension [0, a] x [0, L2] x ... [0, Ld] and the region II of dimension 
[a, L\\ x [0, Li\ x ... [0, Ld]- We are mainly concerned with the case where L\ — ► 00 
which implies that region II is opened. The case where L\ remains finite will be 
discussed in a later section. The two and three dimensional rectangular pistons are 
illustrated in Figure 1. Casimir force acting on the piston due to massless scalar field 
with periodic (P), Dirichlct (D) and Neumann (N) boundary conditions (b.c), as 
well as electromagnetic field with perfect electric conductor (PEC) b.c. and perfect 
magnetic conductor (PMC) b.c. will be computed. 

We first recall some basic definitions. The finite temperature Casimir energy 
is naively the sum of the zero temperature Casimir energy (jl.lj) and the thermal 
correction: 

E Cas = E° Cas + A£ Cas = E° Cas +T £ log (l - e— ' T ) = -Tlog J] 1 _ /T 
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Figure 1: The two and three dimensional rectangular pistons 



For the piston system, the Casimir energy is the sum of the Casimir energies of 
regions I and II and the Casimir energy of the region outside the rectangular cavity: 

(2- 1 ) ^Cas = ^Cas + ^Cas + ^Cas ■ 

Although each of these terms is divergent, we are going to see below that the 
Casimir force acting on the piston defined by 

( 2 - 2 ) ^Cas = --g^E Cas = (^Cas + ^Cas) 

is finite without any renormalization. The term Eq^ s has drop out from the last 
equality of (|2.2[) since it is independent of a. 

Notice that E l Cas = E Cas (a, L 2) ..., L d ) and Egks = e Cob(Li - a, L 2 , ■ ■ ■ , L d ), 
where Eq^Li, L 2 , ■ ■ . , L d ) is the Casimir energy of the rectangular cavity [0, L{\ x 
[0, L 2 ]x. . . [0, L d }. As was proved in [5S1[3D], the Casimir energy Eca. B (Li, L 2 , ■ ■ ■ , L d ) 
for massless scalar field with Pb.c, Db.c. and Nb.c. and for electromagnetic field 
with PEC b.c. and PMC b.c. are related to each other by the following linear 
relations: 

d 

(2.3) E%? (L l7 ...,L d )=2- d Y / (Tlf- ] £ E^(2L mi , . . . , 2L mj ), 

j=l l<mi<...<mj<cZ 

(2.4) 

d 

EcSs?{Li, L d ) = (d - l)E® aa (Li, . . . , L d ) + E E^ s (Li, . . . , Lj-i,Lj + i, ...,L d ) 

3=1 

and 

(2.5) E™ C (L U . . . , L d ) =J2U ~ 1) E ^ (*5 L mi ,...,L mj ). 

j — 2 l<m\<...<rrij<d 

Notice that when the dimension d is equal to one, there are no electromagnetic field 
satisfying either PEC b.c. or PMC b.c; whereas the Casimir energy for massless 
scalar field with Db.c. and Nb.c. is related to the Casimir energy for massless 
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scalar field with Pb.c. by the simple relation (a) — (l/2)EQ as (2a). In the 

following discussions, when we mention massless scalar field with Db.c. and Nb.c. or 
electromagnetic field with PEC b.c. and PMC b.c, we implicitly assume that the 
space dimension d is larger than two. 

For a massless scalar field with Pb.c. in a rectangular cavity [0, Li] x [0, L 2 ] x 
... [0, L d ], the eigenfrequencies are 



E, 2nk 1 



{ki,...,k d ) e Z d . 



Let 



(2.6) ^f s (A;L 1 ,...,L (i )=i £ 



E, 2nk 7 



\U 



keZ d \{0} 



exp —A 




,L d ) 



be the A-dependent zero temperature Casimir energy. Here the function K(\\ L±, . . . . L d ) 
is defined by 



K{X;L U 
Using the formula 



we find that 
(2.7) 



Z exp 

kez d \{o} 



/ 2nkj x 2 



— / r(w)z""'dw, Rc c> 0, 

ZlTl J c^ioQ 



-, r c+ioo I d r„ , 



n 2 



= 5-7 / r(«j)A Z d — 



d 

Re c> -, 



where Zd(s; ai, . . . , a d ) is the homogeneous Epstein zeta function |31[ 132] defined 
by the infinite series 



Z d (s;ai,...,a d ) = I X^'*^ 

kez d \{o} \j=i 



(2.8) 



when Re s > d/2. The function Z d ( ^; . . . , |^ has a meromorphic continua- 
tion to C with a simple pole at to = d with residue 
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By shifting the contour of integration in (|2.7j) from the line Re w = c, c > d/2 to a 
line Re w = — 2 — e, e > 0, we find that 

(2.9) 

d 



r(d + i) 



3=1 



1 



2 d 7T d/2 T (d) 

Using the functional equation (see e.g. [5T1 152"! 1551 154"! 155] ) 
(2.10) K- s T{s)Z d {s;a u ...,a d ) 



1 2tt 



2tt 



2' Li 



1 



r I d \ ( 'I 1 

2 / \2 ai a d 



we find that the cut-off dependent finite temperature Casimir energy is given by 
(2.11) 



^Cas(A;£l> ■ ■ ■ >Ld) 



T(d+1) 

2 d 7r d/2 T (|) 



d 



A d 1 + i?cas, rcg(^l! • • ■ > L d) + 0(A), 



where the regularized finite temperature Casimir energy -E<? as reg (£i, . . . , Lg) is 
given by (see [2"9"]V 



Cas, rcg C^l ? ■ * ' ? 



2^ 



^d I o — ! -^l) ' ' * ? c 



(2.12) 



T ^ log 

k6Z d \{0} 



1 



CX P 



2tt 
T 7 



d 



d+i 



0;T, 



-Tlog 



2tt 



T 



1 

' ' id, 

The cut-off dependent finite temperature Casimir energy for massless scalar field 
with Db.c. and Nb.c. and electromagnetic field with PEC b.c. and PMC b.c. can 
be computed using formulas (|2.3|) , (|2.4p and (|2.5p . It is easy to deduce from (|2.11[) , 
(|2.3[) . (|2.4p and (|2.5[) that in all cases, the cut-off dependent Casimir energy can be 
written as a sum of the A — > + divergent term and the regularized term: 

£'Cas(A; L\, . . . , Ld) = Ecas. div(A; Li, . . . , Ld) + Ecas, reg(il 5 • ■ ■ i L d ) + 0(A). 

Although the divergent term for the massless scalar field with Pb.c. only depends 
on the bulk volume, we can see from (|2.3p , (|2.4p and (|2.5p that for massless scalar 
field with Db.c. and Nb.c. and electromagnetic field with PEC b.c. and PMC b.c, 
the divergent terms depend on the area of lower dimensional hypersurfaces. As 
mentioned in the introduction, it has been argued by several authors [3j [H [5j [6l 
[3 [8] that these hypersurface divergences cannot be removed by renormalization of 
physical parameters. However, we notice that regarded as a function of L\, the 
divergent term is linear in L%, i.e. 

(2.13) Ecas, Div(A; L\, . . . , Ld) — ei(A, L 2 , . . . , Ld)L\ + eo(A; L 2 , ■ ■ . , Ld)- 

Here we want to remark that the above calculations is an idealization of the more 
physical description of the interaction between the quantum field and an external 
potential, which goes to zero away from the boundary surfaces. For a more rigorous 
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treatment, one can consider the approach of [5] where the imposed boundary con- 
ditions are approximated by adding a <5-type background potential to the quantum 
field which is concentrated on the boundaries. In our idealization where the piston 
is assumed to be a perfect rigid rectangular cavity partitioned into two regions by a 
perfect rigid piston with negligible thickness, the Casimir force acting on the piston 
would be independent of the approach used to compute the Casimir energy, as is 
already pointed out by [9]. However, we would like to emphasize that we do not 
claim to have renormalized the Casmir energy in a physical way. Instead, our point 
is that for an idealized piston, the divergence terms which plague the calculations 
of Casimir energy can be ignored since they would not contribute to the Casimir 
force acting on the piston (see next section). But this only holds in the idealized 
situation. For more details about physical situations that would invalidate our 
assumptions, one can refer to [5] and a more recent preprint [36] , 



3. Casimir force on ^-dimensional rectangular pistons 

Since the A — * + divergent term of the Casimir energy for a d-dimensional 
rectangular cavity is linear in L\ (|2.13|) . it follows that for the Casimir energy of 
the piston system (|2.1[) . the A — > + divergent term depends only on L\, and not 
on a. Therefore the A — > + divergent part of the Casimir energy do not contribute 
to the Casimir force acting on the piston. Consequently, the Casimir force can be 
computed by using only the regularized Casimir energy: 



-Fbas(a; L 2 , . . . , L d ) = - lim — \ E Cas , ro g (a, . . . , L d ) + E Cas , reg(£i - a, ... , L d ) 



For massless scalar field with Pb.c, we obtain from eq. (|A.3[) that the regularized 
Casimir force acting on the piston due to the vacuum fluctuation of the field in 
region II is given by 



(3.1) 




(3.2) 




and the total Casimir force acting on the piston is 



(3.3) 
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for massless scalar field with Pb.c. For massless scalar field with Db.c. and Nb.c, 
(f2~2"]) and (231) give 

(3.4) 

d 

F c D a f (a; La, ... , L d ) = 2~ d+1 "£(T) d ~ j £ F Ca,(2a; 2L roi , . . . , 2L m ^J. 

j — 1 2<mi<...<nij_i<d 

For j = 1, the second summation is a single term F(? as (2a). Using (|3.3p . we find that 
the Casimir force on the piston for massless scalar field with Db.c. and Nb.c. can 
be written respectively as 



(3.5) 



og / d / , \ 2 



Fg as (a;L 2 ,.,.,L d ) = -*TY: £ £ £ J 1 + ( 2lT ) 



fel = l (fc 2 ,...,fcd)GN ti - 1 i = -oo \j=2 x J 

x exp —2ixk\a 



\ 



d , , N 2 



£ 

3=2 



(2ff) 2 



and 



(3.6) F& s (a;L 2 ,...,L d ) 



T 
2a 



^£ £ 

fci=i (fc 2 ,...,fc d ,0e[(Nu{o})' i - 1 xZ]\{o} 



£^) 2 + (2/T) 2 j expf^fcxa £ 



For electromagnetic field under PEC b.c. and PMC b.c, (|2 .4[) and (|2.5|) give 
(3.7) 

d 

F c£ C ( a 'i Li,..., L d ) = {d - l)F(? as (a; L 2 , . .. ,L d ) +£ F(? as (a; L 2 , . . . , Lj-%,Lj + x, . . ., L d ); 

3=2 

and 
(3.8) 

F C^ C ( a ' L li ■■■i L d) = £(■?' - 1) £ ^cL(a;£roi. • • • i^mj'-i)' 

j=2 2<mi<...<T?ij_i<c! 

Using these formulas and the formula (|3.5[) . we conclude that 
(3.9) 



^PEC/PMC 
Cas 



(a;L 2 



' 7rT £ £ A PBC / PMC (fc 2 ,...,/ 

fei=l (fc 2 ,...,fc d )e(Nu{o})< i - 1 



£ (£(^y + ( 2 ^) 2 ] exp I -2-Kkxa 



Z=-oo \j=2 



\ 



3=2 
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where 



d-1. 



(3.10) A^(k 2 ,...,k d )= {l, 

0, 



if all fcj, i = 2, . . . , d are nonzero, 
if exactly one of the ki is zero, 
if more than one of the ki is zero; 



(3.11) A PMC (k 2l . . • , kd) = j, if exactly j of the ki, i = 2, . . . , d is nonzero; 

and 6 PEC = 1 if and only if d = 2 and 5 PMC = 0. 

Notice that since the Bessel function K u (z) is positive for any positive z, we 
immediately obtain from (|3.3[) , (|3.5[) , (|3.6[) and (13. 9|) that /or massless scalar field 
with periodic, Dirichlet and Neumann boundary conditions and for electromagnetic 
field with PEC and PMC boundary conditions, the Casimir force acting on the 
piston always has negative sign, and therefore is attractive. This holds for any 
dimension and size of the piston as well as for any temperature. 

By taking derivative of the Casimir force with respect to a, one finds from (|3 . 3[) . 
()3.5p . (|3.6|) and (|3.9p that the derivative of the Casimir force with respect to a is 
always positive. Therefore, we can conclude that for either massless scalar field 
with periodic, Dirichlet or Neumann boundary conditions or electromagnetic field 
with PEC and PMC boundary conditions, the Casimir force is always an increasing 
function of a. Since the Casimir force is always negative, this implies that the 
magnitude of the Casimir force is always decreased when a is increased. To the 
best of our knowledge, this is the first time such results are obtained analytically. 

4. Asymptotic behaviors of the Casimir force for large and small 

plate separations 

The formulas (|3. 3[1 . (|3.5p . (|3.6p and (|3.9p are ideal for studying the large-a and 
small Lj, 2 < j < d behaviors of the Casimir force. In particular, the leading order 
terms of the Casimir force when a — ► oo are given respectively by 



Cas 



(a;L 2 



T 

47rTmin{L 2 

a 



,L-\T}e 



— 2-kcl min{L 9 



,L-\T} 



(4.1) 



F cL(a; L 2, ■■■,L d ) ttT, 



^Cas( a ; L 2 , ■ ■ ■ , L d ) 

(4.2) 



T 
~2a 



J. r-l OT^l „-27ramin{L 2 1 L d 1 ,2T} 



7rTmin{L 2 ,....L d ,2T}e 



T 



Fc£ («; Lt,...,L d )~- 5^ - - TrTmin 

la 





d 


f . j 












d 






exp 


—2ira min < 






\\ 




J=2 




A 


1=2 



(4.3) F^ c (a;L 2: ...,L d )^-nTmm{L^,....L^}e- 2 ^ in ^ ' >. 

The case where T = will be considered in the next section. These asymptotic 
behaviors imply that when a — > oo, the magnitude of the Casimir force tends to zero. 
Moreover, for massless scalar field with Pb.c. and Nb.c. and electromagnetic field 
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with PEC b.c. in d = 2 dimension, the Casimir force tends to zero polynomially, in 
the order a -1 . However, for massless scalar field with Db.c. and for electromagnetic 
field with PEC b.c. and PMC b.c, the magnitude of the Casimir force decays to 
zero exponentially fast. 

For the Casimir force when the plate separation a is small, (|A.4|) and (|3. 2|) give 

us 

(4.4) 



F Cas( a 'j 1*2,..., L d ) 



7r 2 a' 



d+1 



j=2 



4-7T 

d + 4 

a 2 



x Kd-2 



d 



E E ^ E(^) s 

fci=i (fe2,...,fed,Oez d \{o} \i=2 



T 



\ 3=2 



T 



When a — > + , the first term tends to infinity, of order a - . The second term is 
0(a) and the third term tends to exponentially fast. In other words, when the 
plate separation a is small, the leading term of the Casimir force for massless scalar 
field with Pb.c. is 



(4.5) F^ s (a;L 2 ,...,L d ) 



^±1 rf+1 

7T 2 i 1 



i=2 



l) + 0(a°). 



For massless scalar field with Db.c. and Nb.c. and electromagnetic field with PEC 
b.c. and PMC b.c, the corresponding expression for the Casimir force when the 
separation a is small can be found by substituting (|4.4[) into (|3.4[) . (|3.7|) and (|3.8[) . 
In particular, we find that when a is small, the asymptotic behaviors of the Casimir 
force are given respectively by 



Cas 



(a; L 2 , . . . , L d ) 



2 d + 



1 d 

Lew 



< ,ir(^i)C*(i + i) 



i±i 1+1 

7r 2 a^ +i 



S,- + 0(a°), 



and 

^PEC/PMC , T T \ 

where 

(4.6) 5,- = 



J=1 7r _ 2-a-' +i 



E 

2<mi<...<mj_i<d 



Z/ mi . . . Z/ m ^_ 1 . 



Therefore for massless scalar field with Db.c. or Nb.c, the small-a leading term of 
the Casimir force is 

d 

- I / .,!// . I i 



aT(*±l)C R (d + l) 



2 d+i 7r ^±V+i 



J=2 
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which is 2 d+1 times smaller than the leading term for the massless scalar field with 
Pb.c. For electromagnetic field with PEC b.c. and PMC b.c, the small-a leading 
term is 



d 
3=2 



which is times smaller than the leading term for the massless scalar field with 
Pb.c. and d — 1 times larger than the leading term for the massless scalar field with 
Db.c. 

Notice that these leading order terms are independent of the temperature T. 
Therefore we conclude that at low temperature, the effect of temperature to Casimir 
force is insignificant when the plate separation a is small. On the other hand, the 
leading order terms also show that the Casimir force tends to — oo when a — > + . 
Combining with the fact that the Casimir force tends to as a — > oo and the 
Casimir force is an increasing function of a, we conclude that for massless scalar 
field with Pb.c, Db.c. and Nb.c. and electromagnetic field with PEC b.c. and PMC 
b.c, the Casimir force always increases from —oo to as the separation a increases 
from to oo. 



5. LOW TEMPERATURE AND HlGH TEMPERATURE EXPANSIONS OF THE CASIMIR 

FORCE 

As discussed in |29j . the low temperature expansion of the regularized Casimir 
energy inside a rectangular cavity is just the zero temperature Casimir energy plus 
the temperature correction. In the case of massless scalar field with Pb.c, it is 
given explicitly by (|2.12|) . Together with l|3.2j) and (|A. 1[) . we find that the low 
temperature expansion of the Casimir force acting on the piston is given by 



(5.1) 



L 2 ,..., L d ) = F^=°(a; L 2 , . . . , L d ) + J £ £ 

i=i k 1= i (kz,...,kd)ez d - 



-1/2 



J=2 



L, 



exp 



2irl 



\ 



ttT 2 



2tj: y: 

i=i (k 2 ,...,k 4 )ez d - 1 \{0} 



\ 



d 'fc x2 



3=2 



d 'k-^ 2 " 



V, 
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for massless scalar field with Pb.c; and by 



(5.2) 



F Cas (a; L a , ... , L d ) = F^°(a; L 2 , . . . , L d ) + J >T E 



Z=l fe 1 = l ( fc2 



E *? 

,fe<i)G(NU{0})< 1 - 1 



x A(fc 2 



.ttT 2 




-1/2 



exp 



^E E A te, 

i=l (/c 2 ,...,fcd)G(Nu{0})' J - 1 \{0} 



2 v J 



#1 



7rZ 

T 



\ 3=2 



for massless scalar field with Db.c. and Nb.c. and electromagnetic held with PEC 
b.c. and PMC b.c. Here 



A D (fc 2 ,...,fe d ) = 1 



if and only if all ki, i = 2, . . . , d are nonzero; 



A N (k 2l . . . , fed) = 1 and A PEC (k 2 , ...,k d ) and K PMC (k 2l ...,k d ) are given by rfXTUD 
and (ETUI) respectively; 6 N = I; S D = S PMC ee and 6 PEC = 1 if and only if 
d = 2. These formulas show that for massless scalar field with Pb.c. and Nb.c, 
the temperature correction is of order 0(T 2 ) when T<1; but for massless scalar 
field with Db.c. and for electromagnetic field with PEC b.c. and PMC b.c, the 
temperature correction terms go to zero exponentially fast when T — ► + . The 
main contribution to the low temperature Casimir force comes from the first term 
of (|5.ip or (|5.2[) , which is the zero temperature contribution. For massless scalar 
field with Pb.c, it is given by 



Cas 



(a; Li, . . . , L d ) — — 



d+1 



d 
J=2 



C«(d + l) 



[(g) 

2^ 



J=2 



Ai-i — ^ — j L 2 , . . . , L d 



4ir 



d + 4 

a 2 



x E E fc^lE^) 2 

fci=i(fe 2) ...,fe (J )ez«i-i\{o} \j=2 



27rfci 



i=2 



The first term is the leading term and is of order 0(a~ d ~ 1 ) when a — > + . We 
have seen in the previous section that for any finite temperature, this is still the 
leading term when a — > + . The second term is O(a ) and the last term goes to 
zero exponentially fast when a — > + . For massless scalar field with Db.c. and 
Nb.c. and electromagnetic field with PEC b.c. and PMC b.c, the leading behavior 
of the zero temperature Casimir force when a — ► + is also the same as the finite 
temperature case. 
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An alternative expression for the zero temperature Casimir force acting on the 
piston which can be used to study the large a behavior is given by 



j-,P t T=0, T j- \ 

^Cas W L 2,---,^dj 



2 

6^2 " a E E k i\ 

fei=i(te,...,fe ti )ez«i-i\{o} \ 



-<*E E E 

fci=i (fc 2 ,...,fc d )ez<*-i\{o} \j=2 



E 

i=2 

2 



K\ 2irkia. 



\ 3=2 



d N 2' 



\ j=2 



rf /, x2' 



Li 



for massless scalar field with Pb.c. For massless scalar field with Db.c. and Nb.c. and 
electromagnetic field with PEC b.c. and PMC b.c, (pll)) . ((XT)) and (fiTgj) give im- 
mediately 



^Cas (a;L 2l . ■ L d ) 



24a 2 



^E E A(fc 2 ,...,fc rf )fci 

fci=i (fc 2 ,...,fcd)e(Nu{o}) d - 1 \{o} 



\ 



d , kS 2 



3=2 



1« 



\ 



Ell 1 



i=2 



-vrE E A(fc 2 , ■ ■ ■ ,k d ) E ( 7^ ) Ub 27rfc ia 



fci=i (fc 2 ,...,fe d )e(Nu{o})' i - 1 \{o} 
Using the fact that 



,J=2 



as z — > oo, 



\ J=2 



we find that when a> 1, the leading order terms of the zero temperature Casimir 
force are given respectively by 

F Ca S ^2, • • • , id) - 4vra » (mm{L 2 , . . . , L d }J e 



•^Cas {a;L 2 ,...,L d 
^Cas {a;L 2 ,...,L d ) 



ME L 7 2 ] exp(-2™, 

3=2 



\E^ 

\ 3=2 



24a 2 



5 (minjL^ 1 , . . . , e -2™min{L.-\...,L d - 1 } ; 











F c as (a;i 2) ...,i rf ) ~ o 24fl2 Tra 2 


H 












3=2/ 



T?PMC( n . T T 

^Cas (a,^2, ■ ■ ■ 





' . j 














x exp 


— 27ra min < 






[') 






\ 


Z=2 





7ra 2 



(min{Z 2 ,...,L d }) 3 e 



-1 V^ p -2iramin{L 2 1 ,...,L d 1 } 
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We notice a considerable difference between the a — ► oo leading behavior of the 
Casimir force in the zero temperature case and finite temperature case. At zero 
temperature, the leading order terms of the Casimir force for massless scalar field 
with Pb.c. and Nb.c. are of order a~ 2 , in contrast to the order a -1 when the 
temperature is nonzero. For massless scalar field with Db.c. and for electromagnetic 
field with PEC b.c. and PMC b.c, the zero temperature Casimir force also tends 
to zero exponentially fast when a — > oo; but it tends to zero faster than at finite 
temperature. 

Now we consider the high temperature behavior of the Casimir force. From 
(|A.5[) . we find that for massless scalar field with Pb.c, the contribution from region 
I to the Casimir force acting on the piston is given by 



(5.3) F P 6 I as (a;L 2 ,...,L d ) = 



r(^i)Cfl(d+i) 



d+i 
7T 2 



J=2 



T 



d+1 



2 da \ ' a ' L 2 ' ' L d 



3=2 



E E 

(=i kez d \{o} 



x (fc lfl ) 2 + J^ikjLj) 2 Ki 2ttIT 



i=2 



\ 



87ra 



.7=2 



oo oo 



^EE E 

1=1 ki = l (fcs,...,/^)^- 1 

d+2 



(ha) 2 + J2(k J L,y 

3=2 

I 2 k{ 



x ( {k ia y + Y J (k L ] y 

J=2 



K d+ 2 2-kIT 



\ 



J=2 



The leading term is the force due to black body radiation and it is independent of 
a. For the contribution from region II to the Casimir force acting on the piston, 
notice that 



• F Cas( a ;^2, ■ ■ ■ ,L d 



r lim F cL( L i- a '' L 2 ,...,L d ) = - lim Fci(a; L 2 , . . . , L d ) 



Therefore, we find that the black body radiation contributions from regions I and 
II cancel with each other, and the Casimir force acting on the piston is given by 

(5.4) 



^Cas( a ; L 2 , ■ ■ ■ , L d ) — -z-^-Z' d I 0, . r 



T 8,(11 1 

2 a^oo da \ a L 2 L d 



T9 z'(o- 1 1 1 

o" 7T d u > _ j "7 ' ■■■i7" 

2 oa \ a L 2 L d 

d 

4 ' ~* 



j=2 



OO oo 



r-EE E 

i=i fe 1= i (fe 2 ,...,fc d )ez d - 1 



(fcia) 2 + ^(fc 3 Lj) 2 ] Ki ( 2vr/T 

j=2 



\ 



J=2 
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87ra 2 



d 
3=2 



oo oo 



3=2 



=1 fei=l (k2,...,k d 


)GZ d - 1 


d + 2 

















(fcia) 2 +£(W 

3=2 



Notice that the high temperature leading term is linear in T, given by the sum 
of the first two terms. It is called the classical term of the Casimir force. The 
remaining (last two) terms decays exponentially when T — > oo. 

Since the sign of the Stefan-Boltzmann term in (|5.3|) is positive, the contribution 
to the Casimir force from region I will become very repulsive at high temperature. 
However, the black body radiation of region II creates a counter force to cancel this 
repulsive force and we find that the leading order term of the Casimir force acting 
on the piston at high temperature is an attractive force (see (I5.5[) ) of order T. 

Now let us compute explicitly the classical term of the Casimir force for massless 
scalar field with Pb.c. Using (|A.2|) . we find that the classical term of the Casimir 
force acting on the piston for massless scalar field with Pb.c. is 



(5.5) 



FP.classical / T 
Gas {a;L 2 ,.. 

\ fcl=l 



■,L d ) 

£ 

(fc 2 ,...,fc d )ez<i-i\{o} 



\ 



d ( h x 2 



3=2 



exp 



-27rfcia 



\ 



d 'k sy ' 



3=2 



A drawback with the formula (|5.5[) is that it does not manifest the small a behavior 
of the classical term. Using (|A.2I) again, we obtain an alternative formula for the 
classical term 



(5.6) 



F, 



P, classical 



Cas 



a;L 2 ,...,L d ) = -T\--^-T[ -)„,•(,/) 



3=2 



a - 4na 2 



d 
3=2 



E E 

fc I =l(fc 3 ,...,fc < ,)eZ'»-l\{0} \3=2 



27rf 



3=2 



Zd-l ( 2' ' ' ' ' Ld 



for massless scalar field with Pb.c, which is suitable for studying small a behavior. 

The high temperature expansion for Casimir force acting on the piston due 
to massless scalar field with Db.c. and Nb.c. and electromagnetic field with PEC 
b.c. and PMC b.c. can be obtained using (|3 .4(1 . (|3.7p and (|3.8j) . For the contribution 
from region I, (|3.4| . (|3.7|) . (|3.8p and (|5.3p show that leading order terms of the 
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Figure 2: Left: The dependence of the Casimir force Fq^ (a; b, c) on the plate separation 
a when T = 0, 1, 5, 10. Right: The deviation of the Casimir force from the classical term 
when a = 0.05, b = 1, c = 1. 



Casimir force from region I are given by 



F, 



D/N,l 



Cas 



(a;L 2 



E 



(T) d - 



.r(^)c«(j + i) 



O(T) 



and 



PEC/PMCS 



Cas 



(T) d - J C2j 



S, 



+i 



0(T) 



respectively, where Sj is the partial hypersurface area defined by (|4.6p . Notice that 
in contrast to the Pb.c. case, now we have temperature corrections of order T l for all 
2 < I < d besides the leading Stefan-Boltzmann term of order T d+1 . However, all 
these terms are independent of a and therefore they cancel with the corresponding 
terms from region II. The high temperature leading order term of the Casimir force 
acting on the piston is still the classical term of order T and is given explicitly by 



(5.7) 



Fclassical/ . j 
Cas \ a i L 2 

A(fc 2 , ...,k d ) 



E 



fci=i (fc 2 ,...,fc d )e(Nu{o}) d - 1 \{o} 




By taking the derivative of (|5.5p and (|5.7p with respect to a, we find that the 
classical term of the Casimir force is also an increasing function of a. In other words, 
the magnitude of the classical term decreases with increasing a. An alternative 
expression suitable for studying the classical term of the Casimir force at small 
plate separation can be derived using (|3.4[) , (|3- T[) , f|3 .8[) and (|5.6|) . 

We would like to remark that the thermal Casimir effect for electromagnetic field 
in a three dimensional rectangular cavity is considered in [40] . In that paper, the 
authors obtain Casimir force that depends linearly in T in the high temperature 
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regime by subtracting the terms 



d 




(5.8) 



£ (T) d - j (2j-d-l) 



so that in the large volume limit, the Casimir force tends to zero. In our approach, 
no subtraction is required since the contribution from the terms (|5.8[) automatically 
cancel. 

Before ending this section, we would like to comment on the nomenclature 'clas- 
sical term' for the high temperature leading term. To restore the constants H, c and 
k into the expressions for the Casimir force, we replace T by kT/(hc) everywhere 
and multiply the overall expression by he. Notice that a term with order T J will 
be accompanied with H J . Since in the high temperature limit, the Casimir force 
acting on the piston is 0(T), this implies that it has a finite classical h — *■ + limit. 
Moreover, the leading term linear in T gives the classical limit and therefore it is 
called the classical term. 

In Appendix B, we compute the explicit expressions for the zero temperature 
Casimir force when d = 1, 2 and 3 from the formulas in this section. The results are 
found to agree with the existing results in [9l [lGj [TJJ [15j [16j HZ] ■ Explicit formulas 
for the low and high temperature expansions of the Casimir force are also given. In 
Figure 2, we show the dependence of the Casimir force Fq^ (a;b,c) on the plate 
separation a and temperature T for three dimensional electromagnetic field with 



Consider the case where the piston is confined in a closed rectangular cavity. 
More precisely, we consider the scenario where the L\ — > oo limit is not taken. In 
showing that the Casimir force is divergence free in the beginning of Section 3, we 
only use the fact that the divergent part of the Casimir energy in a rectangular 
cavity depends linearly on L\ (without passing to the limit L\ — ► oo). Therefore 
in the present case where L\ is finite, it is still true that the Casimir force acting 
on the piston -Fcas (a; L\, Ld) is free of divergence. On the other hand, the 
derivation of (13. 3p from (IA.3j) shows that after some cancelations, the Casimir force 
i 7 Cas(a; Li, . . . , Ld) for finite L\ case can be written in terms of the Casimir force 
F Cas (a; L 2 ,..., L d ) of infinite L\ case by 

(6.1) F Cas (a;Li, ...,L d )= F Cas (a;L 2 , ...,L d )- F Cas (Li - a;L 2 , ■ ■ ■ ,L d ). 

This is true for either massless scalar field with Pb.c, Db.c. or Nb.c. or electro- 
magnetic field with PEC b.c. or PMC b.c. A trivial consequence of (|6.1[) is that 
the Casimir force acting on a piston placed in a closed rectangular cavity is equal 
to zero when the piston is placed exactly in the middle of the cavity. On the other 
hand, since we have shown that in all the cases we consider, Fc &s {a; L 2 , . . . , L d ) is 
an increasing function of a, this implies that if the piston is placed closer to the left 
hand side (i.e. a < L\ — a), then the Casimir force on the piston Fc as (a; L\, . . . , Ld) 
is negative and tends to pull it to the left. In other words, the Casimir force acting 
on a piston which is placed inside a closed rectangular cavity always try to collapse 
the piston to the nearer end. 



PEC b.c. 



6. Comment on piston inside a closed cavity 
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The low and high temperature expansions of the Casimir force i 7 bas(a; L\, . . . , Ld) 
in the case of finite L\ can easily be computed using the formula (|6.1j) and the re- 
sults for the case of infinite L\ in the previous sections. Therefore we omit them 
here. 



7. Casimir force density when the distances of some pairs of 
transversal plates are large 

In this section, we want to study the asymptotic behavior of the Casimir force 
when d—p > pairs of the transversal plates are large, i.e., a, Li, . . . ,L P <C L p+ \ = 
. . . = Ld for some p > 1. More precisely, we are going to derive the limit 

^bas(a; L\,..., Ld) 



Fcas{p, d; a;L 2 ,..., L p ) 



lim 

L„+i = ... = Ld— >oo 



L p +i ■ ■ ■ Ld 

which we call the Casimir force densitjQ. We consider the low temperature expan- 
sion and high temperature expansion separately. The computations are similar as 
in the previous sections. We only write down the final answers. Details can be 
found in the preprint version of this paper posted in the arXiv [H] , 

7.1. Low temperature expansion. In the low temperature regime, the Casimir 
force density acting on the piston is given by 

•^cLfe d; a,L 2 ,..., L d ) = T^ s =0 {p, d; a, L 2 , . . . , L d ) - 2T~^~ 



oc p 

E EE 

(fc 2 ,...,fe p )ezp-i\{0} l=i \j=2 



8ir _ d- P -i 



fcl = l (k 2 ,...,k p )i 





E E E t + E £ 



j=2 



x Kd- 




rf-p+2 
2 



d-p + 2 
7T 2 



I — k\ 



T d - p+2 ( R (d-p + 2) 



for massless scalar field with Pb.c. The expressions of the zero temperature Casimir 
force density which are suitable for studying the small a and large a behaviors are 
given by 



•^Cal °(p,d;a,L 2 , ...,L d ) 



arm 

d+l 
7T 2 



x Z 



E,p-1 



d+l 



Lo 



E^i) 5 
J= 2 




;=2 
i=2 



-i r(^±i) 



27T 2 



J=2 



E E < 

fei=i (fc 2 ,...,fc p )ezp- 1 \{0} 



E( fc ^) 2 ] ; 

3=2 



^When p = 1, this is actually the pressure on the piston when the cross section of the piston 
becomes an infinite hyperplane. 
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and 



f 



P,T=0 
Cas 



(p,d;a,L 2 , ...,L d ) = - 



d-p+ 1 / d-p + 2 



d-p + 2 
7T 2 



2{d-p+l)a 2- 



x K d-p+2 2irkia 



E E 

k 1 = l(k 2 ,...,k p )£ZP- 1 \{0} \j=2 




( R (d-p + 2)a 

d-p + 2 



-d+p-2 



x fc 1 2 Kd-p 2nkia 




— Ana 2 



E E 

fci=i (/c 2) ...,fc p )ezp- 1 \{0} \j=2 




respectively. For massless scalar field with Db.c. and Nb.c. and electromagnetic field 
with PEC b.c. and PMC b.c, the results can be obtained from that for massless 
scalar field with Pb.c. by using 

^c a ^ d ; a; £ 2 ,. 

=2-f +1 E ^Mj + d-p;2a;2L mx ,...,2L mj _ 1 ); 

j — 1 2<?7ii <.. .<rrij_i <p 

(7.2) 

Fcfs C (P, d ; a;L 2 ,..., L p ) = (d- 1)^^, d; a; L 2 , . . . , L p ) 

p 

+ E -^CasiP — 1, d — 1; a; L 2 , . . . , Lj-%,Lj + x, . . . , L p ); 
i=2 



(7.3) F^(p,d;a;L 2 , 



■ ^2(d - p + j - 1) ^2 FczstiJ + d-p;a;L mi ,...,L mj ); 

j — 1 2<mi<...<mj <p 



which can be derived from (f57Sj) . (|377|) and (j3~g)) (see [29 ). We notice that when 
p > 2, the leading order terms of the temperature correction to the Casimir force 
density is of order T d ~ p+2 when T < 1 for massless scalar field with Pb.c. and 
Nb.c. However, for massless scalar field with Db.c. and electromagnetic field with 
PEC b.c. and PMC b.c, the temperature correction to the Casimir force density 
decays to zero exponentially fast when T — > 0. The main contribution to the low 
temperature Casimir force density comes from the zero temperature Casimir force 
density which at small plate separation a, has leading order proportional to — a^ d ^ 1 
that goes to negative infinity. At large a, the leading order term is proportional to 
_ fl -d+p-2 f or mass i ess scalar field with Pb.c. and Nb.c; and decays exponentially 
for massless scalar field with Db.c. and electromagnetic field with PEC b.c. and 
PMC b.c. 

In the particular case where p = 1, i.e. when the piston becomes a pair of infinite 
parallel hyperplanes, we find that the low temperature expansion of the Casimir 
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pressure acting on the piston is given by 

dT (4+1) 

v<L{d;a) = y-^Gi^+i)"-^ 1 



d+1 

TT 2 



r (4+1') 



7T 2 



oo oo 



Ta J 



+ — T 2 EE fc 2 ' 2 ^ 

a 2 fc=u=i 

for massless scalar field with Pb.c. For massless scalar field with Db.c. and Nb.c, 
the Casimir pressure is equal to VQ as (d; 2a), i.e. 

r (<tti\ 



Cas 



(d; a) 



dT 4+1) 

/A N d+i CM" + 1 ) a 
(An) — 



2 

d+1 
7T 2 



7T d-2 x — -V X — > d + 2 d-2 / irlk\ 

2 2 a 2 fe=i i= i 2 v Ja / 



For electromagnetic field with PEC b.c. and PMC b.c, the Casimir pressure is 
(d — 1) times larger than that of massless scalar field with Db.c, i.e. 

,D/N, 



vZ C/PMC {d;a) = {d-l)V^{d;a). 



Notice that the zero temperature Casimir pressure on a pair of infinite parallel 
plates is equal to 



(IT (4+1) 



for massless scalar field with Pb.c. ft is 2 



-d-l 



with Db.c. and Nb.c, i.e. V. 



d/n,t=o 

Cas 



times weaker for massless scalar field 



(d;a) 



-d-lpP,T=0 



Cas 



(d; a) . For electromag- 



netic field with PEC b.c. or PMC b.c, the zero temperature Casimir pressure on a 
pair of infinite parallel plates is 



^PEC/PMC 
Cas 



(d; a) = - 



d(d- 1) 
(2 7 r) i 



f 



( R {d+l)a 



-d-l 



These agree with the well known results. For the temperature correction, we find 
that the leading term is of order T d+1 when T < 1. In particular, the leading 
thermal correction to Casimir force for massless scalar field with Dirichlet boundary 
condition is 



T (4+1) 

TT 2 



in agreement with the result in [42] . 



7.2. High temperature expansion. In the high temperature regime, we obtain 
from (|5.4[) the following expansion for the Casimir force density for massless scalar 
field with Pb.c: 



■^CasCP) d;a;L 2 ,...,L p ) 



(fc ia ) 2 + ^(fc,L, 

J=2 



B p (a;L 2 ,...,L p )T + 4 



J=2 



OO oo 



^EE E 

1=1 fe 1= l (k 2 .,...,k p )e2P- 1 



/2 



— d 






2?r/T A 







(ha) 2 + ^2(kjLj 

J=2 
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87ra 2 



3=2 



oo oo 



i=l fe 1= l (fc 2 ,...,fc p )GZP- 1 



3=2 



d + 2 




) K *P 


2tt/T a 




V \ 



(fc ia )2+^(fc,^)2 
3=2 



In the high temperature limit, the last two terms decays exponentially while the 
leading term is the classical term linear in T with coefficient B p (a\ L2, ■ ■ ■ , L p ) 
which can be computed from (|5.6[) . Using (|A.1[) , one gets 



B p (a; L-2, . . . , L p ) 



(d-1). 



7T2 



3=2 

EM 



E E 



fei=i (fc 2 ,...,fc p )ezp- 1 \{0} v=2 
d 



3=2 



Zp-i 




CR(d) 



3=2 



a + 47ra 2 



E^ l j 



d+l 



27rfci 



The corresponding results for massless scalar field with Db.c. and Nb.c. and elec- 
tromagnetic field with PEC b.c. and PMC b.c. can be obtained using (|7.ip . (|7.2p 
and (|7.3|) . Again, we find that in the high temperature limit, the leading order 
term of the Casimir force density is the classical term of order T. It can be very 
negative when the plate separation a is small. 

Consider the particular case where p = 1. The high temperature expansion of 
the Casimir pressure for the massless scalar field with Pb.c. is given by 

V^a) = (I) Ud)a- d T + 4T^a-if / ±0k-i 



7T2 



1=1 fe=l 



00 00 

x K4 (2wlkTa) - Sita-^T^ E E l^k~^Kd±2 (2TrlkTa) . 

1=1 k=l 

For massless scalar field with Db.c. and Nb.c, it is given by 



(7.4) vT(d;a) = (I) CR (d)a- d T + 2-^T^a-i±±lh- 

2 7r2 V ^ i=ifc=i 



OO OO 

d-8 d-2 d + 4 x — ^ X — ^ rf + 2 d-2 

x Ki (AnlkTa) - 2"~ vra ~ — ' 



1 d + 4 x — ^ X — T d + 2 d-1 

-T— 22221— k~—Kd+* (AnlkTa). 



1=1 fe=i 

For electromagnetic field with PEC b.c. and PMC b.c, it is [d — 1) times of (|7.4[) . 
The leading term in (|T.4|) agrees with those given in [42] . 

8. Conclusion 

We consider finite temperature massless scalar field with Pb.c, Db.c. and Nb.c. and 
electromagnetic field with PEC b.c. and PMC b.c. for any space dimension d. For 
a rectangular piston which can be considered as a one-sided open rectangular cav- 
ity divided into two regions, it is shown that one can obtain a finite unambiguous 
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Casimir force acting on the piston. Different exact expressions of the Casimir force 
which are suitable for studying the small and large plate separation limits and low 
and high temperature limits are derived. It is verified analytically that for all the 
cases we considered, although the regularized Casimir force acting on a wall of a 
rectangular cavity can be attractive or repulsive depending on the relative size of 
the cavity, the Casimir force acting on the piston is always an attractive force. 
Moreover, the magnitude of the Casimir force decreases as the separation distance 
between the piston and the opposite wall increases. Another interesting result ob- 
tained in this paper is that at high temperature, the magnitude of the Casimir 
force is found to grow linearly in temperature T. This is in contrast to the result 
for rectangular cavities, where at high temperature the leading term of the Casimir 
force is the Stefan-Boltzmann term of order T d+1 , an order much larger than T. It 
also shows that the Casimir force has a classical h — > limit. On the other hand, we 
also establish that at low temperature, the effect of temperature to the magnitude 
of the Casimir force is insignificant when the plate separation is small. 

We have derived exact expressions for the Casimir force acting on the piston 
which are suitable for studying low and high temperature, and small and large 
plate separation limits. A more detailed numerical study of the results would be 
considered in a future work. The methods used in this paper can be easily extended 
to other quantum fields as well as pistons with arbitrary cross section. In particular, 
it will be interesting to consider massive field, fermionic field or massless field with 
mixed boundary conditions. The later will be a possible candidate for repulsive 
Casimir force. Another possible direction is to study the Casimir piston made of 
dielectric and magnetic materials. 
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Appendix A. Formulas for Epstein zeta function and Casimir energy 

Here we gather the formulas we need for computing the Epstein zeta function 
dHH]). The Chowla-Selberg formula [3J1 [3j0 [33 [38] for Epstein zeta function says 
that 
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Taking derivative at s — 0, we have 
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From this, we obtain the following alternative expressions for the regularized Casimir 
energy Eq^ s rc g (^i' ■ ■ ■ > ^d) lEHl]) f° r massless scalar field with Pb.c.: 
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Appendix B. One, two and three dimensional pistons 

In this appendix, we present the results we obtained in the previous sections to 
special cases with d = 1,2,3. 

B.l. d = 1. As we mention in Section 2, there are no electromagnetic field in 
dimension d = 1. For massless scalar field with Pb.c, we find from (|5.1|) that the 
low temperature expansion of the Casimir force is equal to 
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where F P ^ r ~°(a) is the zero temperature Casimir force given by 



"6a 2 
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For the high temperature expansion of the Casimir force, eqs. (]5.4[) gives 
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This formula can also be directly derived from (|3.3[) . (|B.2[) shows that the classical 
limit of the Casimir force due to massless scalar field with Pb.c. is given by —T/a. 
For massless scalar field with Db.c. and Nb.c, we have 

^Caf(«)=^Cas(2a). 

It is interesting to note that (jB.ip and (|B.2p give us the identity 
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B.2. d — 2. In d = 2 dimension, the Casimir force acting on a piston for massless 
scalar field with Db.c. and electromagnetic field with PMC b.c. coincide; and the 
Casimir force for massless scalar field with Nb.c. and electromagnetic field with 
PEC b.c. coincide. Moreover, 

(B-3) F^ PMC (a; L 2 ) = \ (F^(2a; 2L 2 ) - F^ s (2a)) ; 

F^ s PEC (a; L a ) = ~ (F<? a8 (2a; 2L 2 ) + F&.(2a)) . 



Denote L 2 by 6. Using (|5.1|) . we find that the low temperature expansion of the 
Casimir force acting on the piston for massless scalar field with Pb.c. is 
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When a is small, the zero temperature Casimir force has a representation 
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When a is large, it has another representation 
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The high temperature expansion of the Casimir force can be computed by (|5.4j) 
which gives 

{j . 7 oooc ^ 000000 
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The first term gives the classical limit of the Casimir force for massless scalar field 
with Pb.c. 

For electromagnetic field with PEC b.c. or PMC b.c, we can use eq. (|B.3[) to 
obtain the low temperature expansion and high temperature expansion from the 
corresponding expansions for massless scalar field with Pb.c. For low temperature, 
we have 
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Here ,J P = 1 and S PMC = 0, and the zero temperature Casimir force is given by 
(B.4) 
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For electromagnetic field with PEC b.c. or equivalently massless scalar field with 
Nb.c, eqs. (|B.4|I and (|B.5|) agree with the corresponding formulas in [16] . For 
electromagnetic field with PMC or equivalently massless scalar field with Db.c, 
they agree with the results of [9J [17] . 
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For high temperature, 
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The first term is the classical limit of the Casimir force for electromagnetic field. 

In the b — > oo limit, we find that the pressure on the infinite piston x\ — a due 
to massless scalar field with Pb.c. is 
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For massless scalar field with Db.c. or Nb.c, Vq^ 1 (a) = 'P ( ? as (2a). 



B.3. d = 3. When d = 3, the Casimir force for the electromagnetic field with PEC 
b.c. and PMC b.c. coincide. Moreover, 
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Setting L 2 = b and L 3 = c, we find from (|5 . 1 [) that in the low temperature limit, 
the Casimir force acting on the piston for massless scalar field with Pb.c. is 
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and for massless scalar field with Db.c. and Nb.c. and electromagnetic field with 
PEC b.c. and PMC b.c. is 
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for all (k 2l k 3 ) e (N U {0}) 2 , A PEC (k 2 ,k 3 ) = 2 if both fc 2 and fc 3 are nonzero, 
A PEC (k 2 ,k 3 ) = 1 if exactly one of the k 2 or fc 3 is zero and A PEC (k 2 ,k 3 ) = if 
&2 = &3 = 0. The zero temperature Casimir force for massless scalar field with 
Pb.c. can be expressed as 

Fg=° (a; b,c) = - ^Z 2 (2; b,c) + ™£ £ £ fc? ((fc 20 ) 2 + (MT^ 

fei = l(fc 2 ,fe 3 )GZ2\{0} 



x exp (-^V / (fc2&) 2 + (M 2 ) 



3 



7T 2 &C 7T 2 C Cg(3) 1 4 V" V" f fc2 V R- f 27Tk 2 k 3C 

30a4 90 63 27r C 2 h ^^2^\ k ) K i{ b 

K 2 = l K 3 = l 



+ ^?E E fc i (( fc2fo ) 2 + (fc 3C ) 2 )" 1/2 cx P f-^vtw + (fc 3C ) 2 

a fci=i (fc 2 ,fc 3 )ez 2 \{o} ^ a 



or 



2 / / \ 2 / / / ; \ 2 

K9 



fci=i(fe 2 ,fe 3 )ez 2 \{0} V \ / \ / \ 



fci=i (fc 2 ,fc 3 )ez 2 \{o} V v 7 v 7 / \ v v 

which arc suitable for studying the small a and large a behaviors of the Casimir 
force. For massless scalar field with Db.c. and Nb.c, the corresponding formulas 
are given respectively by 

(B.7) 

DfN,T=0, , v_ 7T 2 &C , b + C IT 2 C Cfl(3)/1 . 1\ Cfl(3) 1 7T 



*cT " = >; 6, c) = -— ^ ± -—3 Cb(3) - ± 



480a 4 167ra 3 ^ tv ; 1440 6 3 32tt \6 2 c 2 / 32tt c 2 96a 2 
1 f> v^/feVr/ f2nk 2 k 3 c\ \ Tib ^ 2 / 2Trkik 2 b\ 



^2 — 1^3 — 1 V fcl=l fc2=l 

OO 

, . /■) / V-rr/.-, — 



+ SE E fci ((fc 2 6) 2 + (fc 3 c) 2 )- 1/2 exp (-^V(W + (M 
a fei=i (fe 2 ,fc 3 )ez 2 \{o} ^ a 



28 



S.C. LIM AND L.P. TEO 



and 
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Here d) D = and 6 N = 1. Eq. (|B.7|) agrees with the results of [HI [16] . For 
electromagnetic field with PEC b.c. or PMC b.c, the large a expansion for the 
Casimir force is still given by (|B.8|) with S PE = 0; whereas the small a expansion 
is 
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For the high temperature expansion of the Casimir force acting on the piston 
due to massless scalar field with Pb.c, we obtain from (|5.4[) the formula 
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The Bessel functions ^3/2(2) and K§j2(z) can be re-expressed by elementary func- 
tions using the formulas 

2 V 2z \ 2 / 2 V 22; V 2 z 

The high temperature expansion of the Casimir force for massless scalar field with 
Db.c. and Nb.c. and electromagnetic field with PEC b.c. or PMC b.c. can be ob- 
tained using (|B.6|) . The first term in (|B.9[) give the classical term of the Casimir 
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force for massless scalar field with Pb.c. It has an alternative expression given by 
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For massless scalar field with Db.c. and Nb.c, the classical term of the Casimir 
force is 
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For electromagnetic field with PEC b.c. or PMC b.c, the classical term is given by 

fc 2 =lfc 3 =l V 7 fci = l (/c 2 ,/c 3 )6Z2\{0} V 7 J 

In the infinite parallel plates limit, i.e. b = c — > oo, the Casimir pressure on the 
plates due to massless scalar field with Pb.c. is 
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^Cas(«) = - + E E I'k-^K^klTa) - sVStt-I £ £ lh'^ K^nklTa). 
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In particular, in the low temperature limit, we find that the the leading order term 
of the Casimir pressure on a pair of infinitely conducting parallel plates is the zero 
temperature Casimir pressure given by the well-known result of Casimir [1], i.e. 

' Gas W~ 24()a 4- 

The leading term of the temperature correction is 

7T 2 T 4 



45 ' 

agreeing with the result of |42j . In the high temperature limit, the leading term of 
the Casimir pressure is the classical term 



nPEC, classical/ \ (r(3)T 

which agrees with the results of [43], |44j 02] . 
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